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GENERIC  FACT 


INTRODUCTION 

The  Generic  Sonar  Model  is  a computer  program  designed  to  aid 
analysts  predict  the  performance  of  sonar  systems.  Included  in  the 
Generic  Sonar  Model’s  library  of  subroutines  is  a modified  version  of 
the  Past  Asymptotic  Coherent  Transmission  Model^  (FACT).  The  modified 
version,  called  Generic  FACT,  was  found  to  be  accurate  in  the  majority 
of  cases  investigated,  but  not  in  all. 

Three  questionable  predictions  by  Generic  FACT  documented  in  this 
paper  pertain  to  bottom  reflection  coefficients,  false  caustics,  and 
subsurface  channels.  The  first  prediction  should  be  corrected  by  re- 
placing the  bottom  reflection  coefficient  subroutine,  and  the  second  by 
changing  one  line  of  FORTRAN  code.  Correcting  the  third  prediction, 
however,  may  require  a substantial  effort. 

One  of  the  desirable  features  of  the  original  FACT  Model  is  that 
predictions  made  it  one  installation  can  be  reproduced  at  another. 
However,  if  the  user  makes  the  changes  recommended  in  this  report,  he 
will  lose  consistancy.  A possible  solution  to  this  dilemma  would  be  to 
Implanent  the  same  changes  at  all  installations,  but  this  is  presently 
impractical. 


GENERIC  SONAR  MODEL 

The  Generic  Sonar  Model  is  a computer  program  designed  to  provide 
systems  developers  and  technologists  with  a comprehensive  modeling 
capability  for  evaluating  the  performance  of  sonar  systems  and  investi- 
gating the  ocean  environment  in  which  they  operate.  Specific  uses  of 
the  Generic  Sonar  Model  include  the  following: 

a.  Analysis  in  which  the  relative  ; :r..  ' of  sonar  systems 

are  compared 

b.  Optimization  of  individual  sonar  designs 

c.  Support  of  sonar  operational  trainers 

d.  Support  of  engagement  models 

e.  Comparison  of  existing  acoustic  models 
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£.  Development  of  Improved  acoustic  models. 

Only  a discussion  of  "e"  follows  because  the  main  thrust  of  this 
report  Is  to  discuss  the  comparison  of  existing  acoustic  models. 

Comparison  of  existing  acoustic  models  Is  accomplished  through 
modular  design,  adhering  to  a strict  Generic  format,  and  using  over- 
laying techniques.  For  example,  an  option  presently  available  In  many 
propagation  programs  allows  one  to  compare  bottom  reflection  coeffi- 
cient models  by  simply  replacing  the  bottom  reflection  coefficient  sub- 
routine. The  Generic  Sonar  Model  allows  the  user  to  replace  entire 
propagation  loss  programs  as  well.  Thus,  a sonar  analyst  can  choose 
appropriate  building  blocks  from  a library  of  subroutines  in  order  to 
obtain  a best  cost/accuracy  configuration  for  the  particular  study  at 
hand. 


The  Generic  Sonar  Model  will  not  be  ready  for  general  distribution 
until  it  is  documented.  However,  local  usage  at  NTJSC  has  uncovered 
some  disturbing  predictions  by  FACT  that  warrant  dissemination  now. 
These  predictions  were  first  discovered  using  Generic  FACT  and  then 
verified  against  the  standard  version’s  output  to  ensure  that  Generic 
modifications  were  not  a source  of  error. 


BOTTOM  REFLECTION  COEFFICIENT  DISCONTINUITIES 

Most  computer  programs  designed  to  model  acoustic  propagation  in 
the  ocean  are  optimized  to  predict  propagation  loss  versus  position  at 
one  frequency,  whereas  sonar  system  analysts  require  propagation  loss 
versus  frequency  with  position  as  the  parameter.  Errors  are  readily 
found  by  investigating  discontinuities  in  plotted  outputs.  As  a result, 
propagation  loss  models  are  "debugged"  at  a constant  frequency,  while 
analysts  who  are  unfamiliar  with  the  idiosyncrasies  o an  acoustic 
model  inadvertantly  optimize  their  systems  at  frequencies  where  the 
computed  propagation  loss  changes  abruptly.  Obviously,  the  design  of  a 
bottom  bounce  sonar  should  not  be  based  on  a bottom  reflection 
coefficient  model  that  displays  significant  discontinuities  at 
frequencies  within  the  operating  band. 

A useful  feature  of  the  Generic  Sonar  Model  la  the  ability  to 
display  the  output  of  its  basic  subroutines  versus  frequency.  Upon 
applying  this  option  to  the  Generic  version  of  the  FACT  bottom  reflec- 
tion coefficient  subroutine,  one  sees  from  figure  1 that  the  sub- 
routine's output  is  riddled  with  discontinuities.  Many  of  these  dis- 
continuities fall  at  sonar  operational  frequencies. 

It  is,  therefore,  recommended  that  the  FACT  bottom  reflection 
coefficient  model  be  smoothed  by  first  resioving  the  code  that  leads 
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to  nonmonotonic  discontinuities  and,  then,  adding  a linear  interpola- 
tion between  the  remaining  breakpoints.  Alternatively,  one  could  use 
the  RAYMODE  bottom  reflection  coefficient  model  that  is  coded  by 
Yarger.^  A sample  output  is  Illustrated  in  figure  2. 

The  bottom  reflection  coefficient  discontinuities  were  apparent  to 
many  FACT  users  who  proceeded  to  find  the  origin  Independently  and 
without  Cenerlc-llke  models.  However,  the  next  questionable  prediction 
to  be  discussed  could  not  have  been  foimd  so  easily  and  serves  to 
Illustrate  the  power  of  the  Generic  Sonar  Ifodel. 


FALSE  f^AUSTICS 

Besides  Generic  FACT,  a second  propagation  loss  subroutine  called 
the  Multipath  Expansion^  resides  in  the  Generic  library.  The  two 
subroutines  agreed  reasonably  well  in  the  majority  of  cases  compared, 
although  FACT  was  approximately  20  times  faster.  Both  were  exercised 
for  the  sound  profile  In  figure  3.  Agreement  was  satisfactory  except 
when  both  source  and  receiver  were  placed  at  a depth  of  640  ft 
(195.072  m) . (See  figures  4 and  5.)  This  section  is  concerned  with 
differences  within  10  fcyd  (9.144  km)  in  range;,  the  next  section  will 
discuss  the  more  significant  differences  beyond  10  kyd  (9.144  km). 

As  stated  earlier,  abrupt  changes  in  propagation  loss  tend  to 
iiidicate  error.  Originally,  Generic  FACT  did  not  show  the  bump  of  4 
kyd  (3,6576  km),  aiid  so  it  was  thought  that  the  Multipath  Expansion  was 
at  fault.  The  corresponding  ray  diagram  (figure  6),  on  the  other  hand, 
revealed  a atrong  focusing  effect  at  the  point  in  question,  suggesting 
that  the  incr''"''s  in  intensity  was  real. 

The  Generic  Sonar  Model  is  structured  so  that  individual  contribu- 
tions to  the  acoustic  pressure,  eigenrays,  can  be  examined  individ- 
ually. It  WHS  soon  discovered  chat  Generic  FACT  generated  a family  of 
eigenrays  corresponding  to  a false  caustic  not  generated  by  the  Multi- 
path Expansion.  This  family  of  relatively  high  intensity  masked  the 
caustic  phenocenon. 

The  reader  may  wish  Co  omit  the  remainder  of  this  section  which 
requires  a detailed  kjiowledge  of  PACT. 

According  to  the  NUSC  listing  of  FACT,^  the  FOBTRAS  card 

15  IF  ((THETC.LE.THMAX).AND.(rriETC.GE.THMIN))  GO  TO  (30,35),  IGTYP 
IMST0R93 

is  responsible  for  the  questionable  eigenrays. 
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Here 


THETC  is  the  angle  at  which  a caustic  intersects  depth  Y(K2), 

THMIN,  THMAX  are  angxxlar  bounds  of  a group  of  rays,  and 

IGTYP  equals  2,  indicating  that  the  range  of  the  rays  in  the 
group  is  fit  with  the  parabola, 

■R  = A(l)  + A(2)<tDTHC  + A(3)*DTHC**2 

A necessary  condition  for  caustics  to  occur  is 

DR/DDTHC  = A(2)  + 2.0*A(3)*DTHC  = 0 

so  that 


DTHC  = -A(2)/(2.*A(3))  INSTOR96 

and 

THETC  " THMIU  + DTHC**2  INSTOR97 

In  order  for  the  caustic  to  belong  to  the  group,  the  inequality 
0 < DThO  < (1TIMAX-I'HMIK)**0.5 

SEUat  be  satisfied.  But,  instead  suppose  that  DTHC  is  small  and  nega- 
tive. Then  INSTOR  93  would  incorrectly  place  the  false  caustic  in  the 
group.  However,  when  the  additional  test  that  DTHC  be  positive  was 
Inserted  in  Generic  FACT,  the  questionable  eigenrays  disappeared  and 
the  true  caustic  became  visable. 


SUBSURFACE  CHANNELS 

FACT  incorporates  asymptotic  techniques  in  order  to  model  low 
frequency  diffraction  effects  not  treated  by  classical  ray  theory.  As 
this  last  example  will  demonstrate,  this  technique  is  inappropriate 
when  the  source  and  receiver  lie  within  a weak  subsurface  channel. 

Again  consider  the  sound  speed  profile  in  figure  3.  A weak  sub- 
surface channel  lies  below  the  100  ft  (30.48  n)  surface  duct  and  ex- 
tends to  a depth  of  1265  ft  (385.572  o) . Figure  7 is  a ray  diagram  for 
the  640  ft  (195.072  m)  depth  source  and  consists  cf  source  angles 
ranging  from  -20°  to  +20°  in  1/2°  steps.  Bottom  bounce  rays  were 
terminated  at  their  second  reflection  in  the  ray  diagrss  for 


4 


TR  5635 


illustrative  purpose.  One  sees  that  the  0°  tay  is  trapped  in  the  sub- 
surface channel  while  the  +1/2°  rays  escaped. 

By  applying  the  eigenray  printout  option  of  the  Generic  Sonar 
Model  to  investigate  the  45  kyd  (41.148  km)  region,  we  found  that  the 
Generic  FACT  was  dominated  by  subsurface  channel  energy  (figure  4) 
while  the  Multipath  Expansion  was  bottom  bounce  limited  (figure  5) . 
Analytic  investigation  supported  by  FFP**  and  RAYMODE^  predictions  shown 
in  figures  8 and  9,  respectively,  substantiate  the  Multipath  Expansion. 

Possible  Improvement  to  Generic  FACT  predictions  could  be  achieved 
by  smoothing  weak  subsurface  channels  out  of  sound  speed  profiles  or  by 
adding  specialized  logic  similar  to  the  FACT  surface  duct  model. 


SUMMARY 

This  report  cautions  the  reader  when  using  FACT  in  certain 
applications.  Possible  Improvements  are  suggested,  but  there  does  not 
seem  to  be  any  practical  way  of  implementing  these  modifications  uni- 
versally. Perhaps  an  in ter laboratory  committee  such  as  the  Panel  on 
Standard  Sonar  Models  (POSSM)  could  be  tasked  to  distribute  standard 
model  updates  at  least  throughout  the  Navy.  Secondly,  this  report 
describes  briefly  a Generic  Sonar  Model  and  how  it  may  be  used  to 
localize  discrepancies  in  other  acoustic  models.  Detailed  documenta- 
tion is  in  preparation. 
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Figure  1.  FACT  Eottrjm  Loss  Versus  Frequency  at  10 
Grazinij  Angle  for  KGS  Province  A 


Figure  2.  RAYMODE  Bottom  Loss  Versus  Frequency  at  lO' 
Grazing  Angle  for  MGS  Province  4 


Figure  4.  Propagatloa  Loss  Versus  Range  Computed  by 

Generic  FACT 


RANGE  IN  KYD 

Figure  5.  Propagatiou  Loss  Versus  Range  Computed  by  the 

IJultlpath  Expansion 


Figure  6.  Ray  Diagram 


Figure  9.  Propagation  Loss  Vcreua  Eange  Computed 
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